Torque Reversal in Accretion-Powered X-ray Pulsars 



Insu Yi , J. Craig Wheeler 2 , and Ethan T. Vishniac 2 
institute for Advanced Study, Princeton, NJ 08540; yi@sns.ias.edu 
2 Astronomy Department, University of Texas, Austin, Texas 78712; ethan@grendel.as.utexas.edu, 

wheel@astro.as.utexas.edu 

ABSTRACT 

Accretion-powered X-ray pulsars 4U 1626-67, GX 1+4, and OAO 1657-415 have 
recently shown puzzling torque reversals. These reversals are characterized by short 
time scales, on the order of days, nearly identical spin-up and spin-down rates, and 
very small changes in X-ray luminosity. We propose that this phenomenon is the result 
of sudden dynamical changes in the accretion disks triggered by a gradual variation of 
mass accretion rates. These sudden torque reversals may occur at a critical accretion 
rate ~ 10 15 — 10 16 <? s _1 when the system makes a transition from (to) a primarily 
Keplerian flow to (from) a substantially sub-Keplerian, radial advective flow in the 
inner disk. For systems near spin equilibrium, the spin-up torques in the Keplerian 
state are slightly larger than the spin-down torques in the advective state, in agreement 
with observation. The abrupt reversals could be a signature of pulsar systems near 
spin equilibrium with the mass accretion rates modulated on a time scale of a year or 
longer near the critical accretion rate. It is interesting that cataclysmic variables and 
black hole soft X-ray transients change their X-ray emission properties at accretion 
rates similar to the pulsars' critical rate. We speculate that the dynamical change 
in pulsar systems shares a common physical origin with white dwarf and black hole 
accretion disk systems. 

Subject headings: accretion, accretion disks — stars: magnetic fields — stars: pulsars: 
general — X-rays: stars 



1. Introduction 

Detailed spin evolution on long time scales has been made available for several accretion- 
powered X-ray pulsars such as 4U 1626-67, GX 1+4, and OAO 1657-415 (for comprehensive 
reviews and data, see e.g. Chakrabarty et al. 1993, Chakrabarty 1995, and references therein). 
These systems have shown puzzling abrupt spin reversals. Before and after the observed torque 
reversals, their spin-up and spin-down torques were largely steady. It is intriguing that the spin-up 
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and spin-down torques before and after reversals are nearly identical. These reversals are quite 
different from the random torque fluctuations seen in some pulsar systems believed to be fed 
by winds (e.g. Nagase 1989, Anzer & Borner 1995, and references therein). The nearly steady 
torques plausibly indicate the existence of ordered accretion disks. The observations indicate that 
the mass accretion rate is gradually modulated with small amplitudes on a time scale of at least a 
year, which is much longer than the typical reversal time scale. 

In the disk-magnetosphere interaction models of the Ghosh-Lamb type (Ghosh & Lamb 
1979ab, Campbell 1992, Yi 1995, Wang 1995), the magnetic torque is a function of the mass 
accretion rate, M. The sign of the torque is reversed (spin- up/down) as M varies when the disk 
inner edge moves past the equilibrium radius at which the torque vanishes (e.g. Lipunov 1992). 
In this picture, however, the torque variation is expected to be smooth and continuous unless M 
varies discontinuously. Although this behavior may be relevant for some smooth torque reversals, 
the observed sudden reversals appear distinct (Chakrabarty 1995). Given the lack of any plausible 
mechanism for discontinuous change of M, which must be tuned to occur near spin equilibrium, 
it is difficult for the existing magnetized disk models to provide an explanation. The observations 
indicate that the mass accretion rates vary little during transition (cf. Chakrabarty 1995), which 
makes the discontinuous change of M unlikely as an explanation. 

In this Letter, we propose a possible explanation for the observed sudden torque reversal. 
We take the neutron star moment of inertia = 10 45 <7 cm 2 , radius i?* = 10 6 cm, and mass 
M* = 1.4Mq. At a cylindrical radius R from the star, the vertical component of the dipole 
magnetic field B Z (R) = B*(R*/R) 3 where B* is the stellar surface field strength. The spin period 
is or angular velocity f2* = 27r/P*. 



2. Keplerian Disk-Magnetosphere Interaction and Inner Region 

In the conventional disk-magnetosphere interaction model of the Ghosh-Lamb type (Ghosh 
Sz Lamb 1979ab), the magnetic field of an accreting neutron star penetrates a geometrically 
thin accretion disk and exerts a magnetic torque. Except in a narrow region near the radius 
where the disk is disrupted, it is assumed that the accretion disk rotation is Keplerian, 
n K (R) = (GMJR 3 ) 1 / 2 , the radial internal pressure gradient is small, the radial drift velocity 
is small, and the disk thickness becomes negligible (Campbell 1992, Yi 1995, Wang 1995 and 
references therein). In such a model, the <f> component of the induction equation in steady state 
gives the azimuthal component of the field 

where we have assumed that the internal viscosity and the magnetic diffusivity are due to a single 
turbulent process. Here we will assume that 7, defined as the ratio of R to the vertical velocity 
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shear length scale ^^/(dv^/dz^, is ~ 1 (e.g. Campbell 1992, Yi 1995). The parameter a is the 
usual viscosity parameter (e.g. Frank et al. 1992) and we take a = 0.3. Although we adopt 
specific values of 7 and a, the constraints on B* could always be rescaled in such a way that the 
exact individual values of the three parameters are not necessary (e.g. Kenyon et al. 1996). The 
inner edge of the disk, where the disk is magnetically disrupted, is at R = R determined by the 
condition that the magnetic torque exceeds the internal torque in the disk (Campbell 1992, Yi 
1995, Wang 1995) which can be expressed as 
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where N c = (j/a)B 2 R 3 , B c = B Z (R = R c ), and R c = {GM^P 2 /At: 2 ) 1 / 3 is the Keplerian corotation 
radius. Integrating the magnetic torque over the disk, and allowing for the angular momentum 
carried by the gas which crosses the inner edge of the disk, the torque exerted on the star by the 
disk is 

7 l-(8/7)(iVi? c )3/2 
N ~6 N ° l-{R o/ R c f/ 2 (2 - 3) 

where N = M(GM^Ro) 1 / 2 (Campbell 1992, Yi 1995, Wang 1995) and N -► as R Q /R C 
approaches the equilibrium point x eq = 0.915. This estimate is based on the assumption that there 
is a negligible flux of angular momentum from R < R Q . Most of the contribution to the torque in 
this type of model comes from field-disk interaction in a narrow region just outside of R Q . Several 
different phenomenological descriptions of disk-magnetosphere interaction cause little practical 
differences to our conclusions (Wang 1995). 

Fig. 1(a) presents an example of the typical smooth torque reversal expected in the models 
of the Ghosh-Lamb type is shown. This example is based on a set of parameters similar to what 
we adopt below for 4U 1626-67. It is clear that any densely sampled spin evolution would reveal 
a gradual and continuous variation of the torque, a generic prediction of the Ghosh-Lamb type 
model. A sudden torque reversal with nearly constant \N\ is hard to explain unless there exists 
an unknown constraining mechanism or the M variation is discontinuous and fine-tuned. It is 
possible that some observed smooth torque reversals (Chakrabarty 1995) could be due to this type 
of reversal model. 

Under suitable conditions the inner parts of the accretion disk may evolve to an optically 
thin, low density state, for example, inside the disk boundary layer around a white dwarf (e.g. 
Paczynski 1991, Narayan and Popham 1993). When M is sufficiently small this transition may 
occur over a broad range of radii within a disk, producing a state in which radiative losses are 
so inefficient that the disk retains a large fraction of the heat generated from the dissipation of 
orbital energy (the 'advective state', cf. Narayan and Yi 1995). The transition to this state is not 
well understood, and may be affected by various external factors, including X-ray irradiation (e.g. 
Meyer & Meyer- Hofmeister 1990) and coronal heating (Meyer & Meyer-Hofmeister 1994). 



Here we will assume that, whatever the details of the transition, the disk will make the jump 
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to a low density state when it becomes possible for it to do so. What is the critical M C rit below 
which the disk becomes hot and optically thin? Observationally, weakly magnetized cataclysmic 
variables generally show a trend in which the X-ray to optical flux ratio decreases as M increases. 
Above a critical rate of ~ 10 16 <7 s -1 , the X-ray emission becomes extremely weak. This has been 
attributed to the transition of the inner region to an optically thin (X-ray emitting) hot accretion 
disk (Patterson & Raymond 1985, Narayan & Popham 1993) at low M. Interestingly, such a 
critical rate is largely consistent with the critical rate ~ 2 x 10 17 a 2 or ~ 2 x 10 16 g s^ 1 for a ~ 0.3 
based on the recently discussed advection-dominated hot accretion disks (Narayan & Yi 1995). 
For a given accretion rate, there exists a critical radius inside of which the disk makes a transition 
to sub-Keplerian while the disk at radii larger than the critical radius remains Keplerian (Narayan 
& Yi 1995). The exact location of this critical radius is not clearly understood yet. We assume 
that most of the magnetic torque contribution comes from the inner region (e.g. Wang 1995) 
which is inside the critical radius, e.g. R cr it ^> R c > R - The relevance of the critical M cr j t for 
the cataclysmic variables becomes more striking when we consider strongly magnetized neutron 
star systems. The inner edge of the accretion disk around a strongly magnetized neutron star 

lies roughly at R Q ~ 5 x 10 8 (5*/10 12 G) 4/7 (M/10 16 5 /s) ' cm which is close to the typical 
white dwarf radius. Therefore, one may ask what would be the effects of the transition, at 
M cr u ~ 10 16 <7 s~ 1 , to the hot optically thin accretion disk on the spin evolution of the pulsars? 



3. Disk Transition and Torque Reversal 



We propose an explanation for the torque reversal based on such a transition. The optically 
thin, hot accretion disk cannot be geometrically thin or Keplerian once its internal pressure 
~ pc 2 s becomes a significant fraction of its orbital energy. After the transition, the disk thickness 
H ~ c s /Qk and radial drift velocity vr ~ ac^/Qx increase. The rotation of the accretion 
disk 0, becomes sub-Keplerian Q < (e.g. Narayan & Yi 1995). In the case of strongly 
magnetized pulsars, direct X-ray observation of such a disk transition is difficult because most of 
the X-ray luminosity, L x ~ GM*M/R*, comes from the surface of the star and the luminosity 
from the disk, which is truncated well above the stellar surface, is limited to ~ GM*M / R Q <C L x . 
The sub-Keplerian disk, however, may have observable dynamical consequences. Once the 
sub-Keplerian rotation is forced on the magnetosphere, the (sub-Keplerian) corotation radius is 
shifted inward and the position of the disk inner edge with respect to the new corotation radius 
is relocated. (In principle, the inner edge of the disk could end up beyond the new corotation 
radius but this would not account for the observed properties of the torque reversal systems.) As 
a result, the magnetic torque changes and there could be a visible change of spin-up/down torque 
on the disk transition time scale. 

If such a transition does occur at a certain critical rate, the most likely time scales 
are the local thermal time scale t t h ~ (aVLx)^ 1 or the disk viscous-thermal time scale 
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t dis ~ R/(ac s ) ~ (R/HjianK)- 1 ~ 10 3 s for a ~ 0.3, 12 ~ 10 9 cm, and M ~ 10 16 # s" 1 (e.g. 
Frank et al. 1992). This is plausible if the transition is mainly driven by the thermal instability 
of the local optically thin region of the accretion flow accompanied by the disk density change on 
the local viscous-thermal time scale (cf. Meyer Sz Meyer-Hofmeister 1994). In fact, the viscous 
time scale for the thin disk is unrealistically long, since the relevant infall time is the one for 
the hot, and thick disk, which will be only a few orbital times. In any case, these time scales 
indicate that the local disk surface density change could occur on a time scale short enough, much 
less than a day, to make the transition appear almost instantaneous. The long term gradual M 
modulation determines the overall evolutionary trend and possibly affects the residuals seen in 
some observations (Cutler et al. 1986, Chakrabarty 1995). The sudden torque reversal does not 
require any short term (~ day) change of M but only that the condition M ~ M C rit be satisfied 
around the time of reversal. 

In order to model the reversal episode, we assume no dynamic vertical motion of the disk 
gas such as winds or outflows. The dominant effect of the transition is to reset the corotation 
radius and disk truncation radius R Q . We take the temperature of the sub-Keplerian hot disk to 
be a constant fraction £ of the local virial temperature, i.e. c 2 ~ £,R 2 Vi 2 K (e.g. Narayan & Yi 
1995). Then the ratio of the disk thickness to radius is H/R ~ £}l 2 and the radial drift velocity 
is vr ~ —a^R^x- Using the radial component of the momentum equation (e.g. Campbell 
1992), we get the sub-Keplerian rotation frequency Q/tl K ~ (l - 5£/2 - a 2 ^ 2 /2) 1 ^ 2 = A. 
where £ — > corresponds to the usual Keplerian limit. We note that vr/RVLk ~ a£ < 1 and 
H/R ~ < 1. Assuming a constant £ or A, after the transition to the sub-Keplerian rotation 
with Q(R) = AQ,k{R) < ^k(R), the corotation radius becomes R' c = A 2 / 3 R C and the new inner 
disk edge is relocated to R' a determined by 
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where TVq = AM(GM if R' ) 1 ^ 2 . The torque on the star after the transition is 

N> _7 1-(%/7){RijR' c ?l 2 
N' 6 1-{R'JR' C ?/ 2 ' 



(3-2) 



The torque vanishes when R' Q /R' C — > x eq = 0.915 as in the Keplerian rotation (eq. 2-3). In our 
discussions, we take a constant A = 0.2 (e.g. Narayan & Yi 1995). The parameters a = 0.3 and 
7 = 1 are assumed to be constant before and after the transition. For the Keplerian disk, the 
equilibrium spin (N = in eq. (2-3)) period is 

P - = [4 - 9s] UJ ItP*g) llO^J \JAM~J \wWs) ' (3 " 3) 



For fi = AQk < Qk, the equilibrium spin period would become longer by a factor 1/A and the 
system begins to evolve toward the newly determined equilibrium after transition. 
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4. Sudden Torque Reversals 



For the spin evolution calculation we integrate the torque equation 

f=-r [TorqueI (4 - 1) 

where Torque = N or N' depending on the physical state of the inner disk. We take a linear 
increase or decrease of M as an approximation to more complex M variations on longer time 
scales. The transition occurs on a time scale <C P*/ \dP*/dt\ before and after the transition. 
The transition at M cr it = 10 15 — W 1G g/s, which is determined by the fits to the observed spin 
evolution, is taken to be instantaneous (cf. t t h, tdi s )- For each torque reversal event, we adjust 
P*, M, and the accretion rate time scale, M / dM/dt . We consider three X-ray pulsars for which 
abrupt torque reversals have been detected. For a given initial spin period P*, a fit gives a set of 
the above parameters. 

4U 1626-67: 4U1626-67 (P* pa 7.7s) was steadily spun-up on a time scale ~ W 4 yr 
(P*/P 2 = -8.54(7) x 10~ 13 s~ 2 ) during 1979-1989. The Keplerian corotation radius 
R c Pa (GM^P^/Att 2 ) 1 / 3 = 6.5 x W 8 cm. The recent BATSE detection of a sudden torque 
reversal to spin-down is puzzling due to its very short time scale and the nearly equal spin- 
up/down rates. The steady spin-down torque suggests that there remains a dynamically stable 
(disk) structure after the sudden reversal (For details of observations, see Chakrabarty 1995). In 
Fig. 1(b), the observed torque reversal event is reproduced by P* = 1.2 x W n G, M = 4 x W 15 g/s, 
and dM/dt = -5 x 10 13 g/s/yr which give R a /R c = 0.58, R'JR' C = 0.95, and R' /R = 0.56. The 
derived accretion rate is slightly lower than the previously quoted values (Chakrabarty 1995), 
and accordingly the derived P* is also lower than the previous estimates (e.g. Pravado et al. 
1979, Kii et al. 1986, Chakrabarty 1995). We note, however, that our estimated values (P* and 
M) can always be rescaled by changing a and 7 (Kenyon et al. 1996). The values of a and 7 
based on first principles are not available. The gradual decrease of M is consistent with the 
observed flux decrease (Mavromatakis 1994). The fit naturally achieves the spin-down torque 
which is slightly smaller than the spin-up torque. The fit requires for the transition to occur at 
M cr it = 3.3 x 10 15 g s _1 . The gradual decrease of the mass accretion rate on a time scale ~ 20yrs 
cannot be due to any viscous or thermal processes operating in the inner region (t t h, tdis)- 

OAO 1657-415: OAO 1657-415 has an observed pulse period P* sa 38s. Recent observed 
spin-up/down torques are P*/P 2 pa —7 x 10~ 12 s~ 2 and pa 2 x 10~ 12 s~ 2 respectively (Chakrabarty 
et al. 1993). For the torque reversal in Fig. 1(c), we get P* = 10 12 G, M = 2.0 x W 16 g/s, and 
dM/dt = — 5 x W 16 g/s/yr. The characteristic M modulation time scale is ~ 0.3yr. The critical 
accretion rate M cr it = 1.1 x 10 16 <?/,s which is somewhat higher than the value required for 4U 
1626-67. 



GX 1+4: GX 1+4 has recently shown a sudden transition from spin-down to spin- up 
around the spin period P* pa 122s (Chakrabarty 1995 and references therein). GX1+4 is peculiar 
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in the sense that despite its very short spin time scale ~ 40yr, the spin equilibrium has not 
been reached. It is likely that M fluctuates or oscillates on a time scale <C 4Qyr near spin 
equilibrium. The spin- down rate in the 1980's, iV-P? 

~ 3.7 x 10 12 s 2 , is not far from the 
1970's spin-up rate, -P*/-P 2 ~ —6.0 x 10~ 12 s -2 . The recent 1994 torque reversal from spin-down 
to spin-up lasted for ~ lOOd. This system also showed very similar spin-down and spin-up rates. 
Although there is no significant spectral change in the hard X-ray emission spectra during the 
spin evolution, the flux appears to be increasing as spin-down torque increases (Chakrabarty 
1995), which is in contradiction to the behavior expected in the Ghosh-Lamb type model (eqs. 
(2-2), (2-3)). The fit shown in Fig. 1(d) corresponds to = 3.2 x 10 12 G, M = 5.0 x 10 15 ^/s, 
and dM/dt = 1.0 x 10 ie g/s/yr. Once again, we note that these parameters could be rescaled 
by different choices of a and 7. It is interesting to observe that the M modulation time scale 
M/(M/dt) ~ 5yr is not far from the detected time scale in the coherent variation of pulse 
frequency residual (Cutler et al. 1986, Chakrabarty 1995). The critical mass accretion rate 
M cr u = 6.5 x 10 15 g/s lies between the two values derived above. The gradual decrease of the 
spin- up torque after reversal is not accurately fit in Fig. 1(d) with the linear M variation. This is 
not surprising given the reported unsteady behavior of the spin-down torque before the reversal 
(Chakrabarty 1995). The observed ~ 300d coherent variation in pulse frequency residual is 
intriguing but such a time scale is not far from the estimated orbital time scale ~ yr (Chakrabarty 
1995). 



5. Summary and Discussions 



The proposed transition is most likely to occur at a critical accretion rate 
M cr it ~ 10 15 — 10 16 <7/s. This suggests an interesting connection between the pulsar systems 
and other compact accretion systems such as cataclysmic variables and black hole soft X-ray 
transients. We speculate that the transitions seen in these systems may be due to a common 
physical mechanism, i.e. disk transition to optically thin hot flow. The model indicates that 
the sudden torque reversal could be a signature of a pulsar system near spin-equilibrium with 

M ~ Merit- 

There are some outstanding issues to look into, (i) The origin of the gradual M modulation 
on a time scale ranging from ~ yr to a few decades remains unknown. The modulation by the 
orbital motion on a time scale ~ yr is plausible in GX 1+4 (Chakrabarty 1995) but is unlikely 
in 4U 1626-67 (Rappaport et al. 1977, Joss et al. 1978, Shinoda et al. 1990, Chakrabarty 
1995) and OAO 1657-415 (Chakrabarty et al. 1993). Several binary precession time scales (e.g. 
Thorne et al. 1986) and the mass flow oscillation time scale (induced by X-ray irradiation, Meyer 
& Meyer-Hofmeister 1990) could be relevant for long time scale (> month) modulations. In 
4U 1626-67 the observed optical and X-ray pulsation frequencies are identical, which has been 
attributed to the reprocessing of X-rays by accretion disk (Ilovaisky et al. 1978, Chester 1979). 
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The direct X-ray irradiation of the secondary star (Hameury et al. 1986) and the disk instability 
(e.g. Smak 1984) may also give rise to M modulations in neutron star systems, (ii) It is important 
to quantitatively understand M cr n and R cr it (cf- Narayan & Yi 1995). (hi) Within our model, the 
observed smooth transitions back to spin-up (seen in OAO 1657-415 and GX 1+4) could result 
from the return from the advective to Keplerian flow. The characteristic time scale for such a 
back transition is likely to be the the viscous disk formation time scale. The observed UV-delay 
time scale, on the order of a day, in cataclysmic variables (Livio & Pringle 1992), may be similar 
to the postulated reverse transition time scale, (iv) It remains unexplained that in GX 1+4 the 
X-ray flux increased during the increase of the spin-down torque (Chakrabarty 1995). If X-rays 
come from the shocked polar accretion, as M decreases, the X-ray emission temperature and the 
apparent flux in a fixed X-ray band could decrease due to radiation drag (e.g. Yi & Vishniac 
1994). A geometrically thick, but optically thin, inner disk is more apt to scatter X-ray emission, 
so that a partially obscured polar cap can actually become more conspicuous as M drops. 

This research was supported in part by NSF Grant 95-28110 (JCW), by the SUAM Foundation 
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Figure Caption 

Figure 1: (a) A typical smooth torque reversal event expected in the Ghosh-Lamb type model. 
This illustrative example assumes B* = 10 12 G, P* = 7.68s and M linearly decreasing from 
7.2 x W 16 g/s to 5.0 x 10 16 £r/s over 10 yrs. The torque passes through N = (eq. (2-3)) and the 
continuous variation should show a wide range of torques before and after the torque reversal. 
Upper panel: time variation of torque. Lower panel: spin-up to spin-down transition. Examples 
of (b) 4U 1626-67, (c) OAO 1657-415, and (d) GX 1+4. The solid lines correspond to models 
described in the text and the dashed lines connect data points adopted from Chakrabarty et al. 
(1993) and Chakrabarty (1995). In OAO 1657-415 and 4U 1626-67, the observed spin periods are 
shown only schematically. 
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ERRATUM 

In the Letter "Torque Reversal in Accretion-Powered X-ray Pulsars" by I. Yi, J. C. Wheeler, 
& E. T. Vishniac (ApJ, 481, L51 [1997]), there are errors in the values of the magnetic fields and 
the mass accretion rates. The correct values are as follows and Figure 1 is revised. 4U 1626-67: 
B t = 6 x 10 n G, M = 2.8 x lO 16 ^, dM/dt = -4.5 x 10 14 g/s/yr, R /R c = 0.760, R' /R' c = 0.987, 
R'jRo = 0.444, M crit = 2.2 x W w g/s. OAO 1657-415: B, = 5x 10 12 G, M = 1.2 x 10 17 g/s, 
dM/dt = -1.2 x 10 17 g/s/yr, M crit = 1.0 x W 17 g/s. GX 1+4: = 1.3 x 10 13 G, M = 4 x 10 16 c//s, 
dM/cft = 1.3 x lO 1 ^/ 5 /^, M cnt = 6 x 10 16 5 /s. 

Figure 1-Corrected 

Torque reversal events in three X-ray pulsar systems. Short dashed lines correspond to 
the smooth torque transition (see text) and thick solid lines correspond to the sudden torque 
transition proposed in the present work. The points connected by long dashed lines are observed 
data points adopted from Chakrabarty et al. (1993) and Chakrabarty (1995). In OAO 1657-415 
and 4U 1626-67, the observed spin periods are shown only schematically. Examples of (a) 4U 
1626-67, (b) OAO 1657-415, and (c) GX 1+4. For the smooth transition events, the parameters 
(£*, M, dM/dt) are (10 12 G, 4.2 x I0 16 g/s, - 2 x 10 15 g/s/yr) for 4U 1626-67, (10 13 G, 2 x W 17 g/s, 
- 5 x 10 17 'g/s/yr) for OAO 1657-415, and (4 x 10 13 G, W 15 g/s, 1.8 x 10 17 'g/s/yr) for GX 1+4. 



